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ABSTRACT 

In order to test the hypothesis that Class I protostellar binary stars are a 
product of ejections during the dynamical decay of non-hierarchical multiple sys- 
tems, we combined the results of new adaptive optics (AO) observations of Class 
I protostars with our previously published AO data to investigate whether Class I 
protostars with a widely separated companion (r>200 AU) are more likely to also 
have a close companion (r<200 AU). In total, we observed 47 embedded young 
stellar objects (YSOs) with either the Subaru natural guide star AO system or 
the Keck laser guide star AO system. We found that targets with a widely sepa- 
rated companion within 5,000 AU are not more likely to have a close companion. 
However, targets with another YSO within a projected separation of 25,000 AU 
are much more likely to have a close companion. Most importantly, every tar- 
get with a close companion has another YSO within a projected separation of 
25,000 AU. We came to the same conclusions after considering a restricted sam- 
ple of targets within 500 pc and close companions wider than 50 AU to minimize 
incompleteness effects. The Orion star forming region was found to have an ex- 
cess of both close binaries and YSOs within 25,000 AU compared to other star 
forming regions. We interpret these observations as strong evidence that many 
close Class I binary stars form via ejections and that many of the ejected stars 
become unbound during the Class I phase. 

Subject headings: instrumentation: adaptive optics, stars: formation, binaries: 
visual 
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Introduction 



It is well known that the majority of main-sequence stars are in binary systems. For 
solar- type stars, the binary fraction o ver all separations is 61% after correcting for incom- 
pleteness (jDuquennoy &: Mayorl Il99ll ). Observations of young binary stars are therefore 
critical to understanding the problem of star formation. To this end, there have been nu- 



merous studies of T Tauri binary properties, summarized by lMathieu et al.l (120001 ) . Overall, 
the T Taur i binary fraction is roughly twice the solar-type main-sequence binary fraction 
( lGhezll200ll ). T Tauri stars also have more bi nary stars at wider separations than solar-type 
main-sequence stars (Kohler fc Leinertlll998l ). The optical survey of 238 T Tauri stars by 
Reipurth fc Zinneckerl (119931 ) found 37 binaries with separations from 1" to 12", but only 1 
triple. They propose that many of these binaries have a companion too close to resolve, and 
that these may be resolved with higher angular resolution observations. 

Surveys of the Class I binary fraction have on ly re cently been conduct ed. Observations 
by IConnellev et al.l (j2008ah . iDuchene et all (12004 ). and Irlaisch et all (120041 ) show that Class 
I YSOs also have a binary fraction excess over solar-type main-sequence stars and have 
a similar binary fracti on and binary s epara tion distribution beyond ~100 AU as T Tauri 
stars. Most recently, IDuchene et al.l (120071 ) reported on a search for binary companions 
among 45 embedded young stars with AO, and found a bin ary frequency of 47%±8% over 
the separation range of 14 AU to 1400 AU. iReipurthl (120001 ) presented a scenario by which 
most stars are born in higher order multiple systems that dynam ically decay. In o rder to 
explain the evolution of the Class I binary separation distribution, IConnellev et al.l (l2008bl ) 
adopted this hypothesis and proposed that the envelope that remains after the Class phase 
helps to keep ejected companions from becoming unbound until the envelope dissipates. 
These scenarios predict that Class I YSOs with widely separated binary companions should 
also have a much closer companion. 

Adaptive optics is a maturing technology that has be en shown to be pa rticularly use- 
ful for investigating the binary properties of young stars. iBeck et al.l ( 120031 ) observed the 
NGC 2024 embedded cluster with AO at K-band. Of 73 stars observed, they found only 3 
binaries and 1 triple, which is consistent with the solar-type main-sequence binary fraction 
within the separation range (145 AU to 950 AU) in which they were sensitiv e to companions. 



Twent y-four nearby pre-main-sequence stars were observed with AO by iBrandeker et al. 



( 120031 ). who found tha t 10 o f these are close binaries and 2 are triple star systems. More 
recently, ICorreia et al.l (120061 ) used AO to observe known binary T Tauri stars, searching for 
higher order multiples. Over the separation range of 12 AU to 2000 AU, they found that 
27% of the total number of systems they observed are higher order multiples, including 7 
triple and 7 quadruple star systems. Laser guide star AO is a new technology that allows 
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AO observations of many more embedded sources than was previously possible with natural 
guide star AO. 

The goal of this study is to explore the possibility that man y widely separated Class I 



binary stars are the result of an ejection. IConnelley et al.l (l2008bl ) predicted that protostars 



with widely separated companions within a projected separation of 5,000 AU should be more 
likely to have a close companion than protostars without a widely separated companion. If 
so, then Class I objects with a closely separated companion may also be more likely to have 
a widely separated companion than Class I objects without a closely separated companion. 
We also consider nearby YSOs that may have been ejected but are now too widely separated 
to be a gravitationally bound companion. In this paper we describe our sample selection 
and observations (Section 2), analyze our results (Section 3) and discuss our findings and 
what we have learned about star formation (Section 4). 



2. Sample Selection & Observations 



This paper uses the AO observations presented in IConnelley et al.l (j2008af ) taken with 
the Subaru Telescope as well as new observations taken with the Keck II telescope. The 
sample of Class I YSOs tha t we ob served with AO was selected from the protostars listed in 
Table 1 of IConnelley et al.l (j2008al ). The Subaru natural guide star adaptive optics (NGS- 
AO) system requires that a R<16 guide star be within ~30" of the target whereas the Keck 
II laser guide star adaptive optics (LGS-AO) system requires that a R<18 tip-tilt refere nce 
star be within ~40" of the target. We used the USNO B1.0 catalog jMonet et alJbooah to 
find visible guide stars within these radii from the list of Class I YSOs for the appropriate 
observing run. We also visually inspected POSS plates to verify that the catalog entry is 
not a false detection or a reflection nebula. In the case of our Keck II observations, we 
also used the Keck Observatory's AO planning tool to ensure that the guide star was in 
a location where the wavefront sensor could use it. In addition to selecting targets based 
on the availability of natural guide stars, we chose sources that are redder than H— K= 1 
to ensure that the source is deeply embedded. We also avoided sources with significant 
reflection nebulosity at K-band since structure in the nebula could interfere with detecting 
a companion. 

We used the 36 element curvature NGS-AO system and CIAO (Murakawa et al. 2004) 
with the 0'.'0213 pixel -1 plate scale on the Subaru telescope on the nights of Nov. 15 and 16, 
2005. All targets were observed at L' (3.8 fxm). The median FWHM angular resolution of 
our Subaru AO L' data is 0'.'138 with a standard deviation of 0'.'035 (the diffraction limit is 
0'.'096). This observing run focused on targets in the Taurus and Orion star forming regions. 
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We were not able to observe all of the targets we intended to observe since a full moon in 
Taurus prevented us from observing many targets in that region and our third night was too 
cloudy to use the AO system. We used the Keck II LGS-AO system and NIRC2 with the 
O'.'Ol pixel -1 plate scale camera during two half-nights, on June 10 and July 9, 2007. Our 
observations were conducted primarily at L'. The median FWHM angular resolution of our 
Keck L' data is 0'.'087 with a standard deviation of 0'.'014, which compares favorably to the 
diffraction limit of ft'078. Many of the close binaries found with our L' observations were 
also observed at K-band. 

Dithering was used for all observations in order to remove the detector flat field effects. 
In the case of L' observations, coadds were used to increase the effective integration time 
per image to ~20 s to increase observing efficiency. For the Subaru observations, three 
images were taken at each dither position to reduce the effect of the overhead involved with 
dithering. For the Keck observations, we stayed at each dither position long enough for the 
low bandwidth wavefront sensor to take a few exposures (which varied with the brightness 
of the tip-tilt reference star). Standard stars that have been observed by UKIRT through 
the MKO filter set (Simons & Tokunaga 2002, Tokunaga & Simons 2002) were selected from 
the UKIRT faint standard star list, and were observed for photometric calibration. 

All data were reduced using the same procedure. A dark frame was made by averaging 
together 10 individual dark frames of the same exposure time as the science data. This dark 
was then subtracted from each target frame. To make the "sky" frame, each dark subtracted 
frame was scaled to have the same median value, then averaged together using a min-max 
rejection. The resulting sky frame was then normalized using the median value of the pixel 
counts. Each dark subtracted (but not scaled) target frame was divided by this normalized 
"sky flat" . The median sky value for each frame was subtracted from each frame to set the 
average background counts in each frame to 0. The images were then aligned and averaged 
together using an average sigma clipping rejection. Since the AO optics are warm, nearly 
all of the L' "sky" brightness is from the telescope and the AO system. The procedure we 
used did not make a true L' flat. However, the U "flats" we used were effective in removing 
the flat-field properties due to the detector, which dominates over the signal of the target in 
the raw data. Figure 1 shows the S'.'OxS'.'O area around each binary star observed with Keck. 



AO images taken with Subaru are presented in IConnelley et al.l (j2008al ) 



We subtracted the PSF of the primary star to look for close and faint companions. 
Nearby field stars made the best PSF reference stars. However, given the narrow field of view 
and the scarcity of bright stars at L', these were rare. In the case of the Subaru observations, 
a radial average of the primary star's PSF was also used as an effective PSF model. This 
technique could not be used for the Keck observations due to the hexagonal pupil, so we 



- 5 - 



often used stars that were observed just before or after the target. However, constructing 
a PSF from other target stars was difficult due to how the AO correction depended on the 
brightness of the tip-tilt reference star and the rotation of the pupil from target to target. 
The Keck PSF was also not point symmetric, so rotating the PSF 180° then subtracting 
resulted in large residuals. 



2.1. Separation and Position Angle Uncertainties 

Since both Keck II and Subaru are altitude-azimuth mounted telescopes, the instrument 
rotator (in the case of Subaru) or field of view rotator (for Keck II) keeps the orientation 
of the field of view fixed and also sets the position angle of the field of view. We have 
noted discrepancies of 15° and 10°, respectively, in the position angles of binaries between 
our Keck and Subaru data of IRAS 21004+7811 and between our Keck and UKIRT data of 
IRAS 18383+0059. Since UKIRT is equatorially mounted, the orientation of the field of view 
is fixed, so the PA measured from the UKIRT image is believed to be accurate. Since these 
were the only cases among our Keck data where the field of view has North left and East 
down, it is likely that we neglected to set the rotator PA before we started the observation. 
All PA measurements from our Keck and Subaru data are only as accurate as the image or 
instrument rotator, which are 10" for Subaru and conservatively estimated as ~36" for Keck 
(private communications). The uncertainties on both separation and PA are based on a 1/4 
pixel centroiding accuracy. 



2.2. Results 



Our observations are summarized in Table 1, including details of the binary companions 
observed in the course of our AO observations. The observations in November 2005 were 
with the Subaru telescope, whereas the June and July 2007 observations were with the 
Keck II telescope. The new binary companions discovered with the Keck II observations are 
IRAS 21004+7811 Ba and Bb, IRAS 21007+4951 B, and IRAS 21445+5712 B (see note in 



Append ix A). Most of the other binary compani ons in Table 1 were found by lConnelley et al. 



(j2008al ). FS Tau Aa and Ab were discov ered bvlSimon et al.l (119921 ). and the l'/82 companion 
to IRAS 06382+1071 was discovered bv lPiche et all Jl995h . 
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3. Analysis 

3.1. Questions to Test the Hypothesis 

To address the goals of this study as stated in Section 1, we asked these questions: 1) 
Are Class I protostars with a widely separated companion more likely to have a close com- 
panion (r<200 AU) than a protostar without a widely separated companion? 2) Are Class 
I protostars with a close companion more likely to have a wide companion? 3) Are Class 
I protostars without a close companion less likely to have an apparently widely separated 
companion? 4) Is there a regional dependence on the fraction of close or wide binary com- 
panions? To answer these questions, we used an outer separation limit of 5,000 AU (within 
which a companion might be gravitationally bound) and 25,000 AU (at such a separation 
an object is unlikely to be gravitationally bound but is close enough that it may have been 
ejected, as described in Section 3.3). The answers to these questions and their implications 
are discussed in section 4. 



3.2. Biases and Target Identification 



The analysis in this paper only includes targets observed with AO. Of the targets we 
observed with AO, 5 are in Taurus, Perseus, Auriga, and Ophi u chus, 1 4 are in Orion, and 
29 are elsewhere. In comparison, the study by IConnelley et al.l ( ]2008al ) includes 40 targets 
in Taurus, Perseus, Auriga, and Ophiuchus, 33 in Orion, and 152 are elsewhere. The results 
of this study are therefore biased in favor of targets in the Orion star forming region and 
against targets in the Taurus, Perseus, Auriga, and Ophiuchus clouds. 

We compared the spectral in dex distributions (Figure 2) of the targets in this study with 
those in IConnelley et al.l (j2008al ) to determine if the AO observed sample of Class I YSOs 
is significantly differen t than their sample of C lass I YSOs. The median spectral index for 
all targets observed by IConnelley et al.l (j2008al ) is +0.79, with a standard deviation of 0.74, 
whereas for this study the median spectral index is +0.50 with a standard deviation of 0.48. 
While this study is thus slightly biased towards more evolved YSOs, the difference is not 
statistically significant. The 2-Sample Kolmogorov-Smirnov test shows that the two spectral 
index distributions are consistent with each other within 3 a confidence limits. The AO 
observed sample is likely to be biased towards more evolved sources since many of the NGS- 
AO guide stars were cases where enough visible light from the embedded YSOs themselves 
was available to guide on. Additional observations of more highly embedded sources are 
needed to address this bias. 
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Not all of the targets observed for this study are Class I YSOs (i.e. have an envelope, are 
not optically visible, etc.), as seen by the fact that some of the targets in this study were used 
as the optical guide star for the AO system. Since IRAS 22376+7455 (#4 and #5)0, IRAS 
05384-0807 (#6), and IRAS 19247+2238 #3 are optically visible stars and are not deeply 
embedded, and thus are not included in this analysis. However, based on their proximity 
to the embedded protostars and optically visible reflection nebulosity, they appear to be 
associated with the dark cloud hosting the Class I objects. These optically visible stars may 
be as young a s the nearby emb edded protostars and may have become optically visible after 
being ejected. iReipurthl (120001 ) proposed that ejections can rapidly change the classification 
of a YSO, and the large number of hierarchical binaries in these three systems (9 in all) 
suggests that ejections played an important role in the dynamical histories of these IRAS 
sources. 



3.3. Binary Separation Ranges and Contamination 

For the purpose of this analysis, a companion with a projected separation less than 
200 AU is considered to be "close", and if greater than 200 AU is considered to be widely 
separated. The key result of this study, that every embedded protostar with a close com- 
panion also has another YSO within a projected separation of 25,000 AU (see below), would 
not be changed if the definition of a "close" companion was as great as 500 AU. 

Ejected stars may be gravitationally bound or unbound 0, and the decline in the num- 
ber of binary companions with separations greater than 1000 AU versus spectral index sug- 
gests that many widely separated companions become unbound during the Class I phase 



(jConnelley et al.ll2008bl ). To demonstrate that this scenario is plausible, consider an object 
ejected at the escape velocity of a 2 M central (binary) star. The ejected object would reach 
a distance of 25,000 AU in about 2 x 10 5 years, which is consistent with the expected Class 
I life time. Objects with a projected separation beyond 5,000 AU are most likely unbound 
(whether or not they were ejected). If these nearby unbound objects are at such distances 
on account of having been ejected, then these objects were likely to have been ejected at a 
velocity higher than the escape velocity, and thus can get to a separation of 25,000 AU in 
less than 2 x 10 5 years. If our hypothesis is correct, then we can expect to find ejected YSOs 



1 The numbers refer to the identifications of the sources in Figure 6 of Connelley et al. (2008a) 

2 An ejection is the dynamical process that changes a non-hierarchical system into a bound or unbound 
hierarchical system. An ejection may or may not lead to the escape of the ejected member from the system. 
A bound hierarchical system is not always stable and may lead to a later ejection and possible escape. 
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(Class I and otherwise) at separations as great as 25,000 AU from our Class I targets. YSOs 
at separations beyond 5,000 AU are considered as being "nearby" but are not considered to 
be a gravitationally bound binary companion. 

The distance to each target (when known ) is given in Tab l e 1. W e used the distance 



estimate to the Orion star forming region by ISandstrom et al.l (120071 ). who measured the 
parallax of a star in the Orion Nebula Cluster to derive a distance estimate of 389^1 pc 
To date, the most commonly used distance to the Orion star forming region was ~470 pc 
(see discussion in Muench et al. 2008). Using either distance estimate has no affect on the 
outcome of our analysis. 

Contamination must always be carefully considered when looking for binary compan- 
ions. In the case of our close binaries, the chance of contamination is negligibly small within 
a projected separation of 200 AU (0.02% for targets in Taurus, Perseus, Auriga, and Orion). 
This value is based on the density of background stars within AL'=4 magnitudes of the 
primary stars in these parts of the sky. For our more widely separated companion stars, 
contamination is a much greater problem since we are also considering YSOs within a pro- 
jected radius of 25,000 AU. We were often able to use the H, K, and U photometry in 



Connelley et al.l (j2008al ) to determine if a candidate YSO within a projected separation of 



25,000 AU has colors consistent with a YSO or a reddened backgrou nd star. In some cases 



we also used the presence of an optical or near-IR reflection nebula (jConnelley et al.l 120071 ) 
to ensure that the candidate YSO is in the same cloud. 

It is important to note that this study is different from many previous studies of binary 
stars that attempt to determine how many companions each star has. We were not trying 
to determine how many companions a star has out to a separation of 25,000 AU, in which 
case we would have had to determine whether each object within that projected separation 
is a YSO or a background star. Rather, we were only looking to see if any object within 
this projected separation has properties (reflection nebula or near-IR colors) consistent with 
a YSO. Thus, contamination by background stars is much less of a problem for this study, 
allowing us to extend the outer search radius to 25,000 AU. 



3.4. Restricted Sample 

Surveys for binary stars often include a correction to compensate for an incomplete 
sensitivity to close and faint binary companions. This is particularly important when data 
were taken under different seeing conditions or targets are at different distances, as is the 
case here. We did not use an incompleteness correction for this analysis since the present 
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study is not investigating the number of binaries in a given separation range, but rather we 
are interested in whether individual targets have a binary companion and, if so, at what 
separation. Thus, in addition to considering all companions to all targets, we considered a 
restricted sample consisting of only targets within 500 pc and discarding binary companions 
with a projected separation of less than 50 AU to minimize incompleteness the effect of our 
incomplete sensitivity to close companions of more distant targets. We chose these values 
since we are typically sensitive to all companions as closely separated as O'.'l (50 AU at 500 pc) 
with a contrast less than AL=3 with our AO data. However, we searched for companions as 
faint as AL=4, where we could typically detect companions only as close as 0'.'34 (50 AU in 
Taurus or Ophiuchus, 160 AU at Orion). Thus, we are typically not sensitive to any "close" 
companions if they are as distant as Orion and are more than AL=3 magnitudes fainter 
than their primary star, but we could detect them in the closer star forming regions. There 
are also two targets in the restricted sample (IRAS 22376+7455 #2 and #6) where our L' 
observations were not deep enough to detect a companion 4 magnitudes fainter than the 
primary star. 



3.5. Triple Systems 

Among the 6 triple systems that we observed with AO, 4 have a companion to the 
primary star and 2 have a companion to the secondar}^. The two cases of a binary sec- 
ondary are IRAS 18383+0059 and IRAS 21004+7811. In the case of IRAS 18383+0059, 
the secondary is just slightly fainter than the primary at L' and is redder, so it is likely to 
be more luminous and thus may actually be the primary star in the system. In the case 
of IRAS 21004+7811, all three stars have similar colors, and the primary is much brighter 
than the companions. Our finding that 4 (or 5) of 6 triples have a binary primary star is 
consistent with the expectation that the least massive star is the most likely to be ejected 
when a non-hierarchical triple system decays via an ejection. 



4. Discussion 

The following discussion of the answers to our four questions uses results based on the 
restricted sample unless otherwise noted. The data supporting this discussion, from both 
the restricted and unrestricted samples, are presented in Table 2. 



3 The primary star is defined as the star that is brightest at L' 
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Question #1: Are Class I protostars with a widely separated companion more likely 
to have a close companion than a protostar without a widely separated companion? Our 
hypothesis predicts that protostars with a widely separated companion should be more likely 
to have a close companion. We find that a nearly equal fraction of Class I protostars, ~ 30%, 
have a close companion whether or not they have a widely separated companion within 
5000 AU. This suggests that the answer to this question is no, contrary to our hypothesis. 
However, when the outer separation limit is increased to 25,000 AU, we found that among 
the protostars in either the whole sample or the restricted sample that do not have a widely 
separated YSO, none have a close companion. In comparison, ~ 47% of those protostars in 
the restricted sample with such a widely separated compani on also have a clo s e com panion. 



Using the binomial confidence estimator in Appendix B.2 of iBrandeker et al.l (120061 ). this is 
a higher percentage with a 99.6% (2.63 a) confidence. The results based on the unrestricted 
sample are consistent with the results from the restricted sample. 

Question j^2: Are Class I protostars with a close companion more likely to have a 
widely separated companion? If our hypothesis is correct, then Class I protostars with a 
close companion should be more likely to have a widely separated companion. Considering 
companions out to a separation of 5000 AU, we find that this is true, but not by a statistically 
significant margin. Only a third of those Class I protostars with a close companion also have a 
companion within 5000 AU. However, we found that every Class I protostar in either sample 
that has a close companion also has another YSO within 25,000 AU. When widely separated 
YSOs are considered as far as 25,000 AU, Class I protostars with a close companion are 
more likely to have a widely separated YSO nearby with a confidence of ~ 3.9 a. Recall that 
companions within 5,000 AU are considered bound, and companions within 25,000 AU to be 
unbound. If our hypothesis is correct, then this result suggests that most (about two-thirds) 
of ejected companions become unbound by the end of the Class I phase. 

Question #3: Are Class I protostars apparently without a close companion less likely to 
have an apparently widely separated YSO? If widely separated companions are created via 
ejections, then those Class I protostars without a close companion should tend to not have a 
widely separated companion. Considering widely separated companions as far as 5,000 AU, 
we find that this is indeed the case in both samples. However, in both samples, we found 
that Class I protostars without a close companion are equally likely to have or not to have 
another YSO within 25,000 AU. In the cases where a Class I protostar does not have a 
resolved close companion but does have another YSO within 25,000 AU, it is possible that 
the Class I protostar does have a very close companion that remains unresolved, or the other 
YSO is within a projected separation of 25,000 AU by chance. 

Question #4 : I s there a regional dependence on the fraction of close or wide binary 
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companions? IConnelley et al.l (l2008bl ) considered the Class I binary separation distribution 
for targets in three groups. Group I includes targets in the Orion star forming region. Group 
II includes targets in the Taurus, Auriga, Perseus, and Ophiuchus clouds. Group III includes 
all targets in other clouds. They found that the Class I binary separation distributions are 
significantly different between Group I and II. Since we observed only a few targets in Group 
II with AO, and since t h e samp le observed in this study is much smaller than the sample 
used by lConnelley et al.l (l2008bl ). we divided the AO observed sample between those targets 
in the Orion star forming region and those that are not in Orion. We stress that the targets 
in Orion are not in the Orion Nebula Cluster, but are distributed throughout the Orion star 
forming region. 

We found that the Orion and non-Orion sub-samples have the same fraction of widely 
separated companions within 5000 AU. However, the Orion sub-sample has a significant 
excess of YSOs within 25,000 AU of the Class I target. We find that there is a 92.5% chance 
(1.4 a) that the Orion and non-Orion sub-samples are different based on these data. We 
found that the Orion sub-sample also has a significant excess of close binary stars versus 
the non-Orion sub-sample (50.0%^'g^ and 15.0%^g 2 ^T°, respectively). The Orion and non- 
Orion sub-samples have a 95.1% chance (1.7 a) of being different based on this data. The 
results from the unrestricted sample are consistent with these findings. The fact that the 
Orion sub-group has these excesses of close binaries and nearby YSOs is consistent with the 
finding that the Orion star forming region (Group I) has a higher Class I binary frequency 
than Groups II or III (17.6%±8.2% vs. 6.1%±3.4% vs. 8.2%±3.4% respectively, Connelley 
et al. 2008b). This finding is also consistent with the hypothesis that young binary stars 
form via ejections from non-hierarchical multiple groups and that the nearby YSOs were 
ejected during the event that created the close binary. It is unclear what cloud properties 
cause the Orion sub-sample to have a higher binary frequency. 



These results are consistent with the scenario presented by iReipurthl (120001 ) in which a 
large fraction of Class protostars form in non-hierarchical multiple systems that dynami- 
cally decay during the Class phase. However, it is important to note that not all widely 
separated binaries are necessarily the result of an ejection. It is possible for a companion 
to have formed at a wide separation, or the presumed companion may be a YSO from else- 
where, or it may be a background star along the line of sight. Our interpretation of the 
observations is that each close binary that we found is the remnant of a multiple system that 
ejected a member, and that the YSOs that are found within 25,000 AU are likely to have 
been ejected. We also found that, among protostars with a close companion, a third have a 
companion within 5,000 AU (which is likely to be gravitationally bound) whereas all have 
a YSO within 25,000 AU (which is not likely to be bound if it is beyond 5,000 AU). This 
result may suggest that about a third of the ejected stars are still bound in the Class I phase 
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and the other two-thirds have become unbound. 



Connelley et al.l (l2008bl ) found that the Class I binary fraction at separations between 
1000 AU and 5000 AU steadily declines throughout the Class I phase (using spectral index as 
a proxy for age), presumably because these widely separated companions become unbound 
as the YSOs age. However, this result is for all Class I YSOs, not just those with a close 
resolved binary companion. Under this scenario, we can expect that Class I YSOs with a close 
companion and a widely separated companion within 5,000 AU will be younger (and thus 
have a higher spectral index) than those Class I YSOs with a close companion but without a 
widely separated companion within 5000 AU (which has presumably become unbound) since 
the dissipation of the envelope allows widely separated companions to become unbound as 
the stars evolve. The spectral indices (measured from 12 fim to 100 fim with IRAS fluxes) 
for those Class I YSOs with both a close companion and a widely separated companion 
within 5,000 AU have a median value of 0.62 ± 0.38, whereas the spectral indices for those 
with both a close companion but without a wide companion within 5,000 AU have a median 
value of 1.03 ± 0.43. The former group only includes 6 IRAS sources, and the latter group 
4 IRAS sources. This result, stron gly limited by small num ber statistics, is contrary to our 



expectations based on the result of I Connelley et al.l (j2008bl ). 



Although the angular resolution of our observations were high enough to resolve binary 
companions as close as 15 AU in the nearest star forming regions, the two closest companions 
we found have projected separations of 36 AU (FS Tau A) and 42 AU (IRAS 21004+7811 B). 
There are several reasons why we did not find binary companions as close as our detection 
limit. Class I YSOs in the closest star forming regions (Taurus, A uriga, Perseus, p Ophi- 



uchus) have a lower binary fraction than Class I YSOs in general (jConnelley et al.ll2008bl ) 
and we only observed 5 targets in these regions. Notably, these two closest binaries are a T 
Tauri star (FS Tau A) and a star with near-IR colors consistent with a T Tauri star (IRAS 
21004+7811, Persi et al. 1988), respectively. Among our bona fide Class I YSOs, the closest 
resolved companion has a projected separation of 65 AU from a source in Orion. Since this 
is near our resolution limit for targets at the distance of Orion, where many of our targets 
are located, it is not surprising we did not find closer companions. 



5. Conclusion 

We observed 47 Class I YSOs using the NGS-AO system on the Subaru telescope and 
the LGS-AO system on the Keck II telescope to search for close binary companions to Class 
I protostars with the goal of testing the hypothesis that protostars with a companion closer 
than 200 AU should be more likely to have a widely separated companion than protostars 
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without a close companion. Our observations found 13 companions with a projected sepa- 
ration less than 200 AU. Contrary to our hypothesis, we found that Class I protostars with 
a close companion are equally likely to have, or not to have, a widely separated compan- 
ion within a projected separation of 5,000 AU. We also found that protostars with such a 
widely separated companion are not significantly more likely to have a close companion. 
However, we found that for every protostar with a close companion, there is another YSO 
within a projected separation of 25,000 AU. We believe that these protostellar binaries and 
widely separated YSOs are the remnants of earlier higher order multiple systems that ejected 
the nearby YSO, now observed at separations up to 25,000 AU. As such, the observations 
support our hypothesis. A third of the close protostellar binaries have a widely separated 
(and presumably ejected) companion within a projected separation of 5,000 AU that may 
be gravitationally bound, while the former companions of the other two-thirds now have a 
projected separation greater than 5,000 AU and are likely to be unbound. The targets in 
the Orion star forming region have an excess of close binary companions and an excess of 
targets with another YSO within 25,000 AU. This is consistent with Orion having a greater 
fraction of close binaries than other star forming regions, and the finding that Class I YSOs 
with a close companion all have another YSO within 25,000 AU. 

This research has made use of the SIMBAD database, operated at CDS, Strasbourg, 
France, and NASA's Astrophysics Data System. This publication makes use of data prod- 
ucts from the Two Micron All Sky Survey, which is a joint project of the University of Mas- 
sachusetts and the Infrared Processing and Analysis Center/ California Institute of Technol- 
ogy, funded by the National Aeronautics and Space Administration and the National Science 
Foundation. The authors wish to thank Dr. Joshua Walawender and Dr. Sean Matt for pro- 
viding helpful commentary during the writing of this paper. The authors wish to recognize 
and acknowledge the very significant cultural role and reverence that the summit of Mauna 
Kea has always had within the Hawaiian community. We are most fortunate to have the 
opportunity to conduct observations from this mountain. 



A. Comments on Selected Targets 



IRAS F04189+2650: The close compani on to FS Tau A was d iscovered by lChen et al. 



(119901 ) . Lunar occultation measureme nts by ISimon et al.l (119921 ) found a separation of 



265 mas ±5 mas and a PA=60° ± 5°. Krist et al. 



Jl998h used WFPC/HST to measure 



a separation of 238.7 mas ±5 mas and a PA=84.4° ± 1.5°. This binary was also observed by 



Hartigan &: Kenyonl (120031 ) with STIS/HST, who measured a separation of 242 mas ±14 mas 



and a PA=93.7° ± 2.5°. The motion of FS Tau Ab relative to FS Tau Aa is shown in Figure 
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3, which includes our observation. The path of FS Tau Ab is remarkably line ar. Since the 



proper motion of FS Tau A is to the SSE (PA=167.5) (IDucourant et al.ll2005l ). the motion 
of FS Tau Ab should be to the NNW if it were a background star. However, its motion is 
too slow and in the wrong direction for it to be a background st ar given the proper motion 



of FS Tau A. Also, spectra taken by lHartigan fc Kenyonl (120031 ) show that FS Tau Ab has 



a late type spectrum with Ha emission and thus is unlikely to be a background star. FS 
Tau B, 19'.'75 to the west, is embedded in a reflection nebula and has H-K and K-L' colors 
consistent with a Class I YSO. 

IRAS 05384—0808: Four stars in this group are binaries. The binary star to the 
west (#6) is optically visible and thus is not a deeply embedded YSO. However, it is in a 
reflection nebula, and so appears to be associated with the embedded stars. 

IRAS 06382+1017: This is a hierarchical triple star. There is another highly embed- 
ded source 13'.'5 to the east that shares a common nebula with the triple star. A bright star 
17'.' 5 to the west is an unrelated foreground star based on its near-IR colors. 

IRAS 07025-1204: The southern group appears to be a hierarchical triple. 

IRAS 18383+0059: This is a hierarchical triple where the companion is itself a close 
binary. However, the brighter star of the close binary (star Ba in Figure 1) is redder and 
only slightly fainter than the primary at L', so it may actually be more luminous. Star Bb 
is found to be much redder (K-L'=3.49) than star Ba (K-L'=1.77). 



IRAS 21004+7811: iPersi et al.l (119881 ) classified this source as a T Tauri star based 



on J through L photometry. This is a hierarchical triple where the companion star is itself 
a 0'.'14 binary with nearly equal brightness at L'. 

IRAS 21007+4951: This deeply embedded (H-K>4, K-L'=5.06) object is found to 
have a very close (0'.'17) companion that is much redder (K-L'=7.78) than the primary star. 

IRAS 21445+5712: We found a possible companion (object B in Figure 1) C'69 away 
that is bluer than the primary star (K-L'=1.50 vs. K-L'=2.56), which is not unusual for 
companion stars. At K-band it is small (FWHM=0 / .'09) but elongated, so it is unlikely to 
be a background star. It is possible, though unlikely, that this is a very small clump of 
nebulosity and not a true companion since it is separated from the reflection nebula. 

IRAS 22376+7455: This nearby (300 pc, Hilton & Lahulla 1995) group has four 
binary stars. The close binary to the NW (#4 in Connelley et al. 2008a) is optically visible 
and was used as the AO guide star. Although it is not a deeply embedded YSO, it is in an 
optical reflection nebula associated with the dark cloud, and thus it is likely to be associated 
with the embedded stars. 
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Fig. 1. — S'.'OxS'.'O field of view around each binary observed using the Keck laser guide star 
AO system at L'. 
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Fig. 2. — Spectral index distribution for the whole sample studied by lConnelley et al.l (j2008al ) 
versus the same distribution for targets in this study. The targets in this study are on average 
younger, but not by a statistically significant margin. The 2 Sample K-S test shows that 
these two distributions are consistent with each other within 3 a confidence limits. 
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Fig. 3. — The motion of FS Tau Ab relative to FS Tau Aa. The uncertainties in our position 
are based on an estimated 1/4 pixel centroiding uncertainty in our image and assume that 
the image rotator accurately put north up. The approximate date of the observation is next 
to each data point. The spacing of the three most recent measurements is prop ortional with 



the tim e interval between these measurements. However, the measurement by ISimon et al. 



(119921 ) in 1989 is too far away, which suggests an error in position angle. The large arrow to 
the right is the proper motion vector of FS Tau A on the sky in the 17.8 year span of the 
observations. FS Tau Aa is located at (0,0) on this figure. 
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Table 1. AO Observed Sample 



IRAS # a a b (J2000.0) <5 b (J2000.0) dist (pc) c AK d AL' d r (") PA UT Date 



03301+3111 




03 


33 


12 


84 


+31 


21 


24 


1 


290(2) 












09Jul07 


F04189+2650 


2 


04 


22 


02 


18 


+26 


57 


30 


5 


140(3) 




2 


67 


0.264+0.003 


106.61+0.76 


09Jul07 


04240+2559 




04 


27 


04 


70 


+26 


06 


16 


3 


140(1) 












14Nov05 


04530+5126 




04 


56 


57 


02 


+51 


30 


50 


9 


none 












14Nov05 


05289-0430E 




05 


31 


27 


09 


-04 


27 


59 


4 


417(6) 




2 


64 


0.288+0.007 


344.62+1.48 


15Nov05 


05289-0430W 




05 


31 


25 


98 


-04 


28 


02 


6 


417(6) 












15Nov05 


05302-0537 




05 


32 


41 


65 


-05 


35 


46 


1 


417(6) 







16 


0.673+0.006 


26.44+0.63 


14Nov05 


05327-0457W 


1 


05 


35 


13 


10 


-04 


55 


52 


5 


417(6) 




1 


23 


0.147+0.007 


82.18+2.90 


15Nov05 


05357-0650 




05 


38 


09 


31 


-06 


49 


16 


6 


417(6) 












15Nov05 


05375-0040 


1 


05 


40 


06 


79 


-00 


38 


38 


1 


417(6) 




1 


68 


6.403+0.007 


278.97+0.07 


14Nov05 


05384-0808 


1 


05 


40 


50 


59 


-08 


05 


48 


7 


417(6) 




1 


81 


0.377+0.005 


141.51+1.13 


15Nov05 


05384-0808 


3 


05 


40 


49 


92 


-08 


06 


08 


4 


417(6) 




1 


89 


0.178+0.007 


355.25+2.40 


15Nov05 


05384-0808 


6 


05 


40 


48 


07 


-08 


05 


58 


7 


417(6) 







73 


0.162+0.005 


314.09+2.64 


15Nov05 


05404-0948 


1 


05 


42 


47 


70 


-09 


47 


22 


2 


417(6) 


3.48 


2 


49 


3.594+0.006 


204.92+0.12 


15Nov05 


05404-0948 


1 


05 


42 


47 


70 


-09 


47 


22 


2 


417(6) 


2.57 


2 


36 


0.162+0.005 


135.32+2.63 


15Nov05 


05513-1024 




05 


53 


42 


55 


-10 


24 


00 


7 


417(6) 












15Nov05 


05555-1405N 


1 


05 


57 


49 


46 


-14 


05 


27 


8 


417(6) 







08 


5.797+0.007 


177.52+0.07 


15Nov05 


05555- 1405N 


1 


05 


57 


49 


46 


-14 


05 


27 


8 


417(6) 




1 


55 


0.209+0.006 


115.31+2.04 


15Nov05 


05555-1405S 


4 


05 


57 


49 


18 


-14 


06 


08 





417(6) 












15Nov05 


06297+1021W 




06 


32 


26 


12 


+10 


19 


18 


4 


900(2) 












15Nov05 


06382+1017 


1 


06 


41 


02 


64 


+10 


15 


02 


4 


800(3) 




2 


46 


1.818+0.006 


336.56+0.23 


14Nov05 


06382+1017 


1 


06 


41 


02 


64 


+10 


15 


02 


4 


800(3) 




1 


85 


0.169+0.007 


17.20+3.12 


14Nov05 


07025-1204 


2 


07 


04 


51 


62 


-12 


09 


29 


9 


1150(7) 







43 


0.345+0.006 


330.32+1.24 


15nov05 


07025-1204 


2 


07 


04 


51 


62 


-12 


09 


29 


9 


1150(7) 




1 


56 


2.364+0.005 


142.76+0.18 


15nov05 


07025-1204 


4 


07 


04 


50 


13 


-12 


09 


15 


7 


1150(7) 




3 


58 


1.488+0.007 


356.69+0.29 


15nov05 


07025-1204 


4 


07 


04 


50 


13 


-12 


09 


15 


7 


1150(7) 




4 


45 


0.638+0.007 


355.13+0.67 


15nov05 


08043-3343 




08 


06 


15 


61 


-33 


52 


19 


5 


1120(4) 












15Nov05 


16235-2416 




16 


26 


34 


17 


-24 


23 


28 


3 


160(1) 












10Jun07 


16316-1540 




16 


34 


29 


31 


-15 


47 


01 


4 


160(1) 












10Jun07 


18250-0351 




18 


27 


39 


53 


-03 


49 


52 





280(8) 












10Jun07 


18275+0040 




18 


30 


06 


17 


+00 


42 


33 


6 


700(9) 












10Jun07 


18383+0059 


3 


18 


40 


51 


73 


+01 


02 


13 


2 


none 


0.81 





05 


2.194+0.003 


280.61+0.56 


10Jun07 


18383+0059 


1 


18 


40 


51 


88 


+01 


02 


12 


8 


none 


4.19 


2 


47 


0.184+0.003 


147.70+1.23 


10Jun07 


19247+2238 


1 


19 


26 


51 


33 


+22 


45 


13 


4 


none 







69 


1.569+0.006 


143.70+0.27 


14Nov05 


19247+2238 


3 


19 


26 


51 


62 


+22 


45 


03 


2 


none 







97 


0.260+0.007 


18.70+1.64 


14Nov05 


19266+0932 
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29 


00 


86 


+09 


38 


42 
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300(3) 
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20355+6343 




20 


36 


22 


86 


+63 


53 
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450(5) 












10Jun07 


20453+6746 
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38 


57 


48 


+57 


09 


37 


6 


500(3) 












14Nov05 


20568+5217 




20 


58 


21 


09 


+52 


29 


27 
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1270(4) 












15Nov05 


21004+7811 


1 


20 


59 


03 


73 


+78 


23 


08 


8 


300(3) 


2.50 


3 


12 


2.515+0.002 


235.35+0.47 


09Jul07 


21004+7811 


2 


20 


59 


03 


13 


+78 


23 


08 


6 


300(3) 


0.35 





01 


0.141+0.003 


119.34+1.50 


09Jul07 


21007+4951 




21 


02 


22 


67 


+50 


03 


08 


1 


700(1) 


3.45 





78 


0.176+0.003 


86.75+1.14 


09Jul07 


21388+5622 


2 


21 


40 


28 


02 


+56 


36 


02 


8 


750(10) 







63 


0.711+0.005 


133.36+0.60 


15Nov05 


21388+5622 


3 


21 


40 


28 


06 


+56 


36 


06 


1 


750(10) 




1 


34 


2.495+0.005 


47.54+0.17 


15Nov05 


21445+5712 




21 


46 


07 


12 


+57 


26 


31 


8 


360(4) 


4.11 


5 


17 


0.694+0.003 


189.68+0.29 


10Jun07 


21454+4718 




21 


47 


20 


66 


+47 


32 


03 


6 


900(1) 












15Nov05 
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Table 1 — Continued 



IRAS 




a b (J2000.0) 


<5 b (J2000.0) 


dist (pc) c 


AK d 


AL' d 


r(") 


PA 


UT Date 


21454+4718 




21 47 20.66 


+47 32 03.6 


900(1) 










09Jul07 


21569+5842 




21 58 35.90 


+58 57 22.8 


250(4) 










09Jul07 


F22324+4024 




22 34 41.01 


+40 40 04.5 


880(3) 










09Jul07 


22376+7455 


1 


22 38 47.02 


+75 11 34.7 


330(1) 




2.22 


0.545+0.007 


175.94+0.78 


15Nov05 


22376+7455 


2 


22 38 43.98 


+75 11 26.9 


330(1) 




0.49 


0.485+0.007 


340.69+0.88 


15Nov05 


22376+7455 


1 


22 38 42.49 


+75 11 45.6 


330(1) 




0.95 


0.695+0.007 


266.78+0.61 


15Nov05 


22376+7455 


6 


22 38 42.83 


+75 11 37.0 


330(1) 




0.16 


1.232+0.005 


142.27+0.35 


15Nov05 


F23591+4748 




00 01 43.25 


+48 05 19.0 


none 


1.64 


1.57 


1.004+0.003 


96.81+0.20 


10Jun07 



a The number of the primary star in the finder charts presented in lConnellev et~ (2008a). 
b Coordinates are from 2MASS, or offset from a field star with 2MASS coordinates 

c Thc estimated distance to each source in parsecs. The citation for the distances estimate is designated by the number 
in the parentheses, and are as follows: 1) Hilton, J., & Lahulla, J. 1995; 2) Educated guess based on proximity to nearby 
objects; 3) Reipurth, A General Catalog of HH Objects, 1999 (http:casa.colorado.edu/hhcat/); 4) Wouterloot & Brand, 1989; 
5) Launhardt & Henning 1997; 6) Menten et al. 2007; 7) Sugitani & Ogura 1995; 8) Bachiller et al. 2001; 9) Zhang et al. 
1988; 10) Battinelli & Capuzzo-Dolcetta 1991 

d The observed magnitude difference between the two binary components in the specified band. 
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Table 2. AO Survey Results, Resticted Sample 







X=5. 


,000 








X=25,000 


Question 


# 
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a 
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% a 


1) Are protostars with a wide companion within X AU more likely to have a close companion 


















than a protostar without a wide companion within X AU? 


















Protostars with a wide companion and nave a close companion: 


3/10 


30 




20 
15 


7% 
.85t 


9/19 


47.3%: 


H3.8% 
-13.4% 


Protostars without a wide companion that have a close companion: 


6/22 


27 


.3%! 


12 
10 


7% 

.2% 


0/13 


0.0%+ 


13.2% 

0.0% 


2) Are protostars with a close companion more likely to have a wide 


















companion within X AU? 










.5% 
.0% 








Protostars with a close companion but without widely separated YSO within X AU 


6/9 


66 


.7% +1 2 7 2 


0/9 


n n%+ 18 - 4% 
0.0/o_ ao% 


Protostars with a close companion and with a widely separated YSO within X AU 
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3) Are protostars without a close companion less likely to have an apparently widely 


















separated YSO within X AU? 


















Protostars without a close companion and without a widely separated YSO within X AU? 
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Protostars without a close companion but with a widely separated YSO within X AU? 
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4) Does the fraction of targets with a wide companion have a regional dependence? 


















Number with a widely separated YSO within X AU in the Orion sub-sample 
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Number with a widely separated YSO within X AU in the non-Orion sub-sample 
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a l <t uncertainties are presented, derived using binomial statistics with the calculator at http://statpages.org/confint.html 



